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In this work we introduce the class of quantum mechanics superpotentials 
W{x) = ge{x)x'^^ and study in details the cases n = and 1. The n = super- 
potential is shown to lead to the known problem of two supersymmetrically related 
Dirac delta potentials (well and barrier). The n = 1 case result in the potentials 
V±{x) = g'^x^ lb 2g|a;|. For V-^ we present the exact ground state solution and study 
the excited states by a variational technic. Starting from the ground state of V- and 
using logarithmic perturbation theory we study the ground states of V+ and also of 
V{x) = g'^x'^ and compare the result got by this new way with other results for this 
last potential in the literature. 



I. INTRODUCTION 

Supersymmetric quantum mechanics (SUSY QM) was first introduced by E. Witten [T] 
[2], as a simplified model (a + 1 dimensional field theory) to study the possibility of SUSY 
breaking. Soon it became a research branch in itself, a way of getting new solutions to 
problems in quantum mechanics [3] |1] |5] |6]. Of particular interest, to our work below, 
we must cite the many papers in the literature [Z] [8] [9] [10] |11] [12] [13] [H] devoted to 
the development of technics for treating the anharmonic oscillator V{x) = bJ^x^ + g'^x'^, and 
other related potentials, which in general do not have exact solutions. 

In this work we present a new simple class of superpotentials in SUSY QM, in the form 
W{x) = ge{x)x'^'"' with n = 0, 1, 2, . . .. The first example of this class, i.e., the case n = 0, 
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was studied long ago in |15j and revisited in |16| and |17| . One of our results is an analytic 
solution for the ground state wave function of the potential V{x) = g'^x^ — 2g\x\, an amazing 
result, considering that analytic solutions do not exist for anharmonic oscillators. Another 
result is a new perturbative solution for the ground state of the potential V{x) = gx^, 
starting from the solution for the potential V{x) = g'^x^ — 2g\x\. Excited states of the 
potentials V±{x) = g'^x'^ ±2g\x\ are also studied by a variational approximation. 

The paper is organized as folows. In Sec. [TT| we make a brief introduction to the well 
known case of superpotentials, which are monomials of odd powers of x, as well as to the 
SUSY breaking ones, which are monomials of even powers of x. More details of these 



solutions can be found in [3] and [H]. In Sec. Ill, we study solutions related to the class 
of simple superpotentials of the form W{x) = ge{x)x'^" {n = 0, 1,2, . . .), where e{x) is the 
sign function. The simple analytical solution for the ground state of the corresponding 
SUSY system is shown, the already known case n = is revised and the case n = 1 is 
studied in more details. The first one is the ilustrative example of the Dirac delta well 
and barrier potentials, which are shown to be SUSY partner potentials associated with 
the superpotential W{x) = ge{x). The second one, W{x) = ge{x)x'^, allows us to find an 



analytical solution for the ground state of the potential V{x) = g x — 2g\x\. In Sec. Ill A 



and Sec. IIIB, we study the excited states of the potentials V±{x) = g'^x^ ± 2g\x\ which are 
derived from W{x) = ge{x)x^. After discussing that exact solutions for the excited states 



cannot be obtained, we apply a variational method (Sec. Ill A) to find approximate solutions 



for the energy levels and the wave functions. In Sec. Ill B , a new perturbative approach to 
the ground state of the potentials V{x) = gx^ and V{x) = g'^x^ + 2g\x\ is presented. Finally, 
a discussion of the results is presented in the Conclusions. 



II. OUR NOTATION AND DEFINITIONS ON SUSY QM 

Let us briefly summarize some main concepts in SUSY QM. For simplicity, we will work 
in a system of units with the Planck's constant set as h = 1 and the particle mass set as 
m = 1/2 (that is 2m = 1). We start by defining the operators A"^ and A: 

A^ = W{x) - ip and A = W{x) + ip , (1) 
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where W{x) is a given function of x and p = —id/dx is the momentum operator. From these 
operators we can construct two hamiltonians: 

H_ = A^A and H+ = AA^ , (2) 

which in terms of p and W{x) resuh in: H± = p^ + V±. The potentials V± are given by the 
equations {W' = dW/dx): 

V±{x) = W{xf ±W'{x) , (3) 

which are Riccati's equations. 

These equations can be understood in two ways. One way is: given W{x), we can define 
the hamiltonians H± with potentials V±{x). The other is: given the potential (or 
V^(x)), by solving one of the Riccati's equations, W{x) can be found, the operators A and 

constructed and the partner potential V^{x) (or V^{x)) can be found. 

The ground state of a SUSY system is defined as the zero energy state of (this is 
a choice; changing the function W{x) — t- —W{x) will change the roles of H and 
As = A^A, its ground state wave function ipoi^) can be got by imposing that it is 
annihilated by the operator A, that is: 

Aij^ix) = (w{x) + ^0 (x) = . 

The solution is given by: 

ij^ix) = Xexp \- ! W{y)dy\ . (4) 



This is a good, physicaly meaningful solution, provided that the function (|4]) is normal- 
izable. Otherwise, a zero energy solution does not exist and SUSY is said to be broken. As 
it is easy to see, superpotentials obeying the rule of being positive (W{x) > 0) for a; > 
and negative (W{x) < 0) for x < shall manifest SUSY. 

Then, starting from W{x), we have two partner hamiltonians, and one of them 
(i7_, in our choice) having a ground state ipQ with energy = and a tower of other 
states: bound states with energies > 0, n = 1, 2, 3, ... or scattering states with energies 
E^ > 0. The hamiltonian has bound energy levels E^_i, n = 1,2,3, .. . with energies 
related to the energies of H_ by the relation: E^_^ = E^ or scattering energies E^ > 0. 



4 



Moreover, the eigenfunctions of and are related according to: 

^ti = {E-y'^' Ai;- (5) 
i^n = {e:.,)-'^' A^^U ■ (6) 

The simplest class of superpotentials manifesting supersymmetry are monomials of odd 
power in x , that is: 

W{x) = gx^"+^ ,n = 0,1,2,... . (7) 
Using the Riccati equation ([s]), we have for the partner potentials: 

V±{x) = W{xf ± W'{x) = ± g^2n + . (8) 



The ground (normalizable) state of H = p"^ + V^, with energy £"0=0 (see Eq. Q) is given 
by: 

The first example of a superpotential of the class ([T]) is W{x) = gx. In this case, the 
associated partner potentials are: 

V±{x) = g^x^ ± g , (10) 

which are simply the potentials of two harmonic oscillators of the same frequency, with a 
constant energy shift g added or subtracted. The ground state of have Eq = 0. Its 
excited states and the states of are given by E~ = E^_^ = 2ng, for n = 1, 2, 3, . . .. 

We will not pursuit the study of this class of superpotentials because they are well known. 
We only mention that the next example of this class, W{x) = gx^, corresponds to the 
potentials V±{x) = g'^x^±3gx'^ and their ground state solution is given by ^ with n = 1. 

On the other hand, the class of superpotentials that are monomials in even powers of 
X, does not give a normalizable zero energy solution to ^ and SUSY is broken. How- 
ever, we can introduce the sign function e{x) and consider superpotentials of the form 
W{x) = ge{x)x'^"'. For this class of superpotentials, a normalizable ground state exists and 
SUSY is not broken. Thus, in the following we study this class of superpotentials, specially 
the n = and n = 1 cases. 
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III. THE CLASS OF SUPERPOTENTIALS OF THE FORM W{x) = ge{x)x'^'^ 

The case n = must be treated separately. So, let us consider the superpotential: 

W{x) = ge{x) , (11) 



where g is a positive constant. For this superpotential (11) the Riccati equations (|3j) give 
the following SUSY partner potentials: 

V±{x) = W{xy ± W'{x) = ±2g6{x) + g"^ , (12) 

where 6{x) is the Dirac delta function. is a delta well, while is a delta barrier, with 
the energy of the ground state displaced by g^. The corresponding Schrodinger equations 
are: 

- il;^"{x) ± 2g6{x)tlj^{x) = [E^ - g^) ^/'^(x) . (13) 

Their solutions are well known [15] [16] |17j . The well iVJ) has a single bound state with 
energy level Eq = 0, binding energy g'^, and wave function given by: 

^o-(x) = v/^e-^l^l . (14) 

All the other eigenstates are plane waves in continuous spectra of energies, the lowest one 
starting with E = g"^. Simple scattering solutions of the well V- and the barrier V"+ can be 
written as: 

^f{x) = A±e'^'' + i3±e-'^" , x < (15) 
^± (x) = C±e^^^ + V^e-'^"" , X > , (16) 



where k = \/ E^ — g^ with E'^ > g^ and the respective constants are related according to 
the bondary conditions ipuiO) = ipiiO) and ip'jjiO) = ip'jiO) ± 2gilj{0) required by the Dirac 
delta potential. 

Summarizing: the hamiltonian if_ has one ground state with energy Eq = and con- 
tinuum of states with energies E^ > g"^ and if+ has a continuum of states with E^ > g"^. 

To see the role of the supersymmetry in this system, let us consider a particle crossing the 
well (or hitting the barrier), coming from x = — oo, such that we can choose V± = 0. With 
the apropriate boundary conditions through x = 0, we can determine B± and C±, getting 
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the scattered and the transmited solutions as functions of the incident amphtudes A±. The 
results can be written as: 

(Tg) 

ijf{x) = A± { e*'^^- + ^ ^e-'^'' } , x<0 (17) 

1-2^) 



k J 

It is easy to verify that the solutions and are related by the supersymmetry 
equations ([S]) and (|6]). For example, by applying the operator A to ^j{x) we get: 



d 
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explicitly showing the manifestation of the supersymmetry of the system. 
Let us now consider the superpotential: 

W{x) = ge{x)x'^ , (19) 

where here also, is a positive constant. The two partner potentials are given by: 

V±{x) = W{xy ± W'{x) = ± 2g\x\ . (20) 

In these potentials a term 6V = ±2gx'^6{x) has been dropped. The reason is that for the 
wave functions involved in this problem its action is null. As the potentials V±{x) — )■ oo for 
X — )■ ±oo, the spectra of H± = p'^ + V± are discrete and their eigenf unctions are normalizable. 
If 6V is treated as a perturbative correction to H±, its action would be non null only if 
f^^dx x'^6{x)\4'{x)\'^ 7^ 0. But this condition requires a wave function that near x = 
behaves like f{x)/x with /(O) 7^ 0, which is non normalizable and is not in the spectra of 
H±. On the other side, treated as part of H±, the term 6V could give non trivial boundary 
conditions for dip/dx ai x = 0. To study this possibility we must integrate the Schrodinger 
equation in the interval x = (— e,e) for e — ?■ 0. A non null effect of 6V only comes if 
J^^dx x'^6{x)ip{x) 7^ 0, which would require a ip{x) behaving like /(x)/x^ with /(O) 7^ 
that is also, out of the spectra of H±. 
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FIG. 1: Partner potentials V-{x) and V+{x) associated with the superpotential = ge{x)x'^. 

A representation of these potentials is given in Fig. [T| As can be seen, V-^ is a single well 
potential and a double well potential symmetric in x. The corresponding Schrodinger 
equations read: 



dx"^ 



+ ± 2g\x\ ) ip^{x) = E^tp^{x) . 



The wave function for the ground state of the double well potential V^{x) = g^x'^ — 2g\x\ 
has energy Eq = and is easily obtained from the equation: 

, s ■, d 

= A^o 



ge[x)x^ + 



The result (already normalized) is given by: 



i)o{x) 



1/3 
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1/6 



(22) 



,2; r(i/3)'/' 

This is an interesting result. As it is well known, exact analytic solutions for the ground 
(or any excited) state of the potentials V{x) = g'^x^ or V{x) = bJ^x^ + g^x'^ cannot be 
obtained. So this exact solution for the potential is somewhat surprising. Another 
characteristic of this solution, is that it represents a single lump centered at x = (which is 
a local maximum of VS) and it is not in the form, as naively expected, of two lumps centered 
at the two symmetric minima, x = ±{l/2gy^^, of V- notwithstanding the fact that, in one 
dimension, any attractive well supports at least a bound state. This happens because the 
"volume" of each well is not big enough to support a bound state (this can be seen in a 
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WKB analysis of the potential, or even more simply, by the Heisenberg uncertaint principle. 
We should only observe that this well size Ax{AEY^'^ is independent of g). 

Let us now look for the excited states solutions. Inspired by the analytic method to solve 



the one-dimensional simple harmonic oscillator and by the form of the solution (22), we try 
a solution of the form 

ij{x) = F(x)e-'''^''/^ . (23) 



Subtituting (23) in the Schrodinger equation (21), it becomes: 



F" - 2ge{x)x'^F'{x) + EF{x) = 



(24) 



For the simple harmonic oscillator, the same steps would lead us to the Hermite equation. 



In our case, we get the equation (24), which is, for a particular choice of parameters, the 
Triconfiuent Heun equation |19|. 



We can go on, look for solutions for the equation (24) through a power series method. 
Assuming that F{x) can be written as: 



(25) 



j=0 



and substituting this expression for F{x) in the differential equation (24), we find: 



3=0 j=0 

Renaming indices and rearranging terms, we have: 



2a2 + Eao + Y t^-^' + 2)0' + 1)«j+2 - '^ge{x){j - l)a,_i + Ea,] = . 
j=i 

Then, given oq e ai, this equation is satisfied if the coefficients aj, j > 2, are given by 
the three terms recursion relations: 

E 



a2 



-ao 



2ge{x){j - 3)aj_3 - Eaj_2 

jU - 1) 



J = 2 
J>3 



(26) 
(27) 



^ In the case of the simple harmonic oscilattor, we supose that the solutions are of the form H{x)e~'' 
and, imposing that those solutions are square integrable, the functions H{x) becomes restricted to be the 
Hermite polynomials ■H„(x^). 
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The corresponding recursion relation for the the harmonic oscilator potential, is a simple 
two terms recursion relation. To get a normalizable solution, we choose the values of E so 
as to terminate the series in a polynomial. In this way we get the set of discretized values of 
the energy spectrum and the corresponding wave functions, that turn up to be the Hermite 
polynomials (see footnote). 



In our case, the recurrence relation (27), is a three terms recurrence relation and there 
is no way of choosing a subset of values of E to terminate the series in polynomials, so as 
to have a normalizable solution. Then, no analytic solution can be found and in the next 



sections we pass to look for approximate solutions. In Sec. |III A]a variational approximation 



is studied and in Sec. [Ill B a perturbative approximation, that will allows us also, to study 
solutions for the potential V{x) = gx'^. 

A. Looking for Approximate Solutions by a Variational Method 

Let us first apply a variational method. The trial function that we are going to use is: 

m 

0(x) = ^«,/,(x), (28) 
i=i 

where j = 1,2, ... ,m. and the coefficients aj G C are the variational parameters. The 
functions fj{x) are chosen to be: 

f^{x) = x^-^e-^'^''/^ . (29) 

This trial function corresponds to the previously used in the power series method, with 
the additional restriction of being a finite polynomial of degree m — 1, instead of an infinite 
series in x. 

For the harmonic oscillator, with a very similar choice of the trial function we would 
find exact solutions. In that case, the variational parameters would be, except for the 
normalization, the coefficients of the Hermite polynomials 'Hn{x'^)- 

Before proceeding, let us consider a convenient change of variables. As can easily be seen, 
by making the rescaling: x — t- g^^^^x it is possible to factor out of the hamiltonians H±, the 
constant g"^^^, that is: 
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So, in the rest of this section, we will work with g = 1 and after finding the energy 
eigenvalues, we can restore the dependence of the energy levels in g by multiplying the 
results by a factor of g"^^^. The restoration of the corresponding wave functions (or trial 
functions), can also be obtained by rescaling x — > g^^^x in the results. 

To go on with the variational method, we construct the expectation value of the energy 
with these trial functions: 



E 



E E«fc«i {fk\H±\fi) 

k=l 1=1 

m m 

E E^fctti {fk\fi) 

k=l 1=1 



(31) 



and minimize E with respect to the parameters ai. This condition gives the system of linear 
equations: 

m 

Y,iiH±)ui - ESui)ai = , (32) 



1=1 



where we used the notation Hki = {fk\H\fi) and Ski = {fk\fi)- The values of E that minimize 
the above system of equations are the eigenvalues of the matrix: 



Mki = {ESu - {H^)ki) 



(33) 



and are obtained by solving the equation det M = 0. The wave functions corresponding to 
each of these eigenvalues are got by substituting the value of E in the linear system above 
and solving for the parameters a^. The matrix elements that we need to construct M^i are: 



kl 



(fklfl) 



dxe sl^l^x''^' ^ 



(34) 



/+00 
dxe-l\^\' [-{1 - m - 2)x^^'-' + 2{l ± l)e{x)x^^'-^] . (35) 
-oo 



0. 



For {k + /) odd, the integrands in (34) and (35) are odd functions and Ski = (-f^±)fc« 
Otherwise, for {k + /) even, we find: 

fc+i-4 

s.H^l ' r(^) (36) 



kl 



-2 



(/- l)(/-2) - (/±l)(A; + /-3) 
(A; + /-3) 



r 



k + i 



(37) 



With these results, the matrix M (33) gets the form: 
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(M±)n 


(M±)3l 



(M±)i3 ... (M±)i„\ 

(M±)22 (M±)24 ... (M±)2„ 

(M±)33 ... (M±)3„^ 



M4 



(38) 



\(M±)^l (M±)^2 (M±)^3 (M±)„4 ... {M±)mm/ 

In this matrix, all elements in positions {k, I), such that (A; + /) is odd are null, while those 



with {k + /) even are given by (33) with Sm and Hm respectively given by (36) and (37). To 
find the energy values we must solve the equation: det M = 0. 

Tables [l] and [ll] show some results found for different number (m) of parameters and for 
g = 1. For different values of g, the values in the Table must be multiplied by a factor of 
g^^^, as observed above. 



TABLE I: Energy values associated with calculated for different numbers of variational pa- 
rameters. 



m 




El 


E2 


^3" 


^4" 


E^ 


Ee 


E- 


1 


0.00000 
















2 


0.00000 


2.04441 














3 


0.00000 


2.04441 


5.76541 












4 


0.00000 


1.97852 


5.76541 


10.00191 










5 


0.00000 


1.97852 


5.54135 


10.00191 


14.94174 








6 


0.00000 


1.97115 


5.54135 


9.49446 


14.94174 


20.37028 






7 


0.00000 


1.97115 


5.51302 


9.49446 


14.06558 


20.37028 


26.29953 




8 


0.00000 


1.96991 


5.51302 


9.41370 


14.06558 


19.02962 


26.29953 


32.64399 


9 


0.00000 


1.96991 


5.50842 


9.41370 


13.90148 


19.02962 


24.43194 


32.64399 


10 


0.00000 


1.96963 


5.50842 


9.39868 


13.90148 


18.73498 


24.43194 


30.18755 



'For this level, the variational method provides the exact solution. 
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TABLE II: Energy values associated with calculated for different numbers of variational pa- 
rameters. 



m 






Et 










Et 


1 


2 


.31447 














2 


2 


.31447 


6.13324 












3 


2 


.04493 


6.13324 


10.54940 










4 


2 


.04493 


5.63655 


10.54940 


15.63469 








r 




i 


.990DD 


0.60600 


9.664 l\) 


-1 r CO Add 

15.63469 


Zl.Zl\)66 






6 


1 


.99066 


5.53888 


9.66470 


14.30956 


21.21933 


27.28556 




7 


1, 


.97666 


5.53888 


9.46567 


14.30956 


19.36916 


27.28556 


33.76558 


8 


1, 


.97666 


5.51611 


9.46567 


13.98107 


19.36916 


24.86727 


33.76558 


9 


1, 


.97235 


5.51611 


9.41524 


13.98107 


18.85787 


24.86727 


30.72924 


10 


1, 


.97235 


5.51007 


9.41524 


13.89369 


18.85787 


24.13659 


30.72924 



The results in Table [T] and |TT] reflect the manifestation of SUSY in the system, at least 
with respect to the equality between the energy levels E~ and E^_^, n > 0, of iJ_ and 
As expected, the ground state energy of is zero and it is not equal to any energy of 
Moreover, for n > 0, increasing the number of variational parameters, we find, mainly for 
the first levels, energies E~ more and more closer to E^_i. 

Therefore, the better the trial we make, the closer we are to satisfy the equality between 
energy levels. Moreover, because the one parameter trial function for the ground state of 
has the same form of the exact (analytical) solution, the value Eq = found is exact 
and the condition of having a zero energy ground state is naturally satisfied. 

Figure |2] shows the first energy levels of H_ and H^. We must remember that the values 
found are better for increasing number of variational parameters and for the lowest levels. 
Thus, for instance, we are supposed to find for the level n = 4 a worse approximation than 
for the level n = 1. 
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FIG. 2: Scheme for the 5 first levels of H- (and 4 first levels of using 6 variational parameters. 

The graphics in Fig. |3] show the approximations for the first levels eigenfunctions of iJ_ 
and respectively. Those approximations were found using 6 variational parameters. 




As expected, we see that the eigenfunctions found have well defined parity, interchanging 
even and odd solutions with even solutions for the ground states. 
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B. Looking for Approximate Solutions by a Logarithmic Perturbation Theory 

We now apply a variant of the logarithmic perturbation theory (LPT) to our problem. 
LPT is explained with more details, for example, in [7], [S], |H] or |11] . 

Starting from the known solution ip^ of VI we can perturbatively obtain the ground state 
of V+, or for example, of the anharmonic potential V{x) = x'^. We start by writting: 

V{x-6) = Vo{x) + 6Vi{x) , (39) 

where: 

Vo{x) = V-{x) = x^ -2\x\ (40) 

Vi{x)=A\x\. (41) 

Observe that V{x; 6 = 1) = and that V{x; 6 = 1/2) = x^. As we only know the ground 
state of V_ we can not go beyond the first order in the Rayleigh-Schrodinger perturbation 
theory. To bypass this difficult we will use the so called logarithmic perturbation theory 
where only the knowledge ipQ is required to calculate the ground state energy level of V{x; 6) 
to any order in 6 (at least numerically). For that aim we consider the perturbed Schrodinger 
equation: 

- ^" + {Vo + 6Vi)^ = (42) 

and write the expansions: 

E = Eo + 6Ei + 6'^E2 + ... (43) 
^ = exp {So + 6Si + 6^82 + ...), (44) 

where Si, S2, etc. are functions and Ei, E2, etc. are numbers to be determined. By 
substituting these expressions in the Schrodinger equation above and equating the terms of 
same powers in 6 we get the set of equations: 

S'^ + S',^ = -Eo + V2 (45) 

S'l + 2S'oS[ = -El + Vi (46) 

S'i + 2S'oS'2 + = -E2 (47) 

SI + 25^5^ + 25^5^ = -E3 (48) 
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Starting with Eq = and 5*0 = — |x|^/3, (that is, \l/o(a;) = ipo — these 
equations can be recursively solved to get Ek and Sk to the desired order in 5. 



Eq. (46) can be rewritten as: 



{S[ exp (2^o))' = (^1 - ^i) exp (2^o) • 



(49) 



By substituting 5*0 = — |xp/3 and Vi = 4|x| in this equation, integrating both sides in the 
interval x = (— oo, +oo) and observing that the integrand of the left side goes exponentially 
to zero at both ends of the integration range, we get for Ei the result: 

(50) 



El 



(^olV^o) 



f^°^ dxe 3 1^1^ 

J — oo 



4 



3V/'r(2/3) 



2.31447 



2 J r(l/3) 

Inserting this result for Ei back into the same equation and integrating now in the interval 
y = (0, x) we get: 



S[{x) = \M^)\-' / dy\My)r[Ei~Vi 

Jo 







2j r(i/3) 



(51) 



9\ ^ 



r(2/3) 
r(i/3) 



r(l/3, 2x73) - r(2/3, 2x^3) 



where Ffa, , are the upper incomplete gamma functions 120]. 



The second order equation (47) can also be written in the form: 



{S'^ exp {2So)y = {-S[^ - E2) exp (25, 



(52) 



Integrating this equation in the interval x = (—00, +00), and observing that the integrand 
of the left side goes to zero at both ends of the integration range, we get E2 as an integral 
over S[: 

(^ol^U^)'l^o) 3 



Eo 



(V'olV'o) 

Substituting (51) in (53), we find 



r(l/3) V3 



dxe-^^''^ S[{xf 



(53) 



Eo 



r(i/3) V3 



2/3 



(Vmi ^in.r (2 2^ .r(2/3) 

\r(l/3)2 ^/''3'3'' +^'/^U'3J ~ r(i/3) '/''3'3''j ' ^ ^ 
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where: 

POO 

Ia{x,y)= I dteH-'^V{x,t)V{y,t) , x>0, y>0, 0<a<l. (55) 



Evaluating the integrals, the expression (54) gives E2 = —0.43817. 



Sumarizing: up to the second order, the ground state energy of V{x] 6) is given by: 

E{6) = Eo + 6Ei + 6^E2 , (56) 

with Eq = 0, Ei = 2.31447 and E2 = -0.43817. 

For (5 = 1, we get for the ground state energy of Vj^, the result: Eq = 1.87630. 

For 6 = 1/2, we find the result Eq^ = 1.04769 for the ground state energy of the the 
quartic anharmonic potential V{x) = x^. This result can be compared with the exact 
one given in |3], noting that our "coupling" constant g is related to their constant g by: 
g2/3 _ ^ly/'^ Thus, multiplying our result by (|)^^^, we find: Eq^ = 0.66000, differing 
of [3] only by about 1.2%. 

On the other hand, comparing the value of E^ found here with the most acurate result of 
the variational method (see Table [ll)) , we see that they differ by about 4.9%, what does not 
seem very good. But, following the suggestion of j8] or |21) . and substituting the expression 



(56) by the corresponding [1, 1] Fade approximant in 6, we find: 



. ^0^, . (57) 

which results (for 5 = 1) in Eq = 1.94605. This result now differs from the result of Table 

I 1 /"^ 

II only by 1.3%. Doing the same for 6=1/2 (and then multiplying by (|) ), we find 
Eq^ = 0.66597, differing from the result of Cooper et all |3] by only 0.03%. A pretty good 
result. 



IV. CONCLUSIONS 

In this paper we studied the class of superpotentials W{x) = e{x)x'^" in SUSY QM. After 
revisiting the case n = we went on studying in details the case W{x) = e{x)x'^. As a 
result we got the exact solution for the ground state of the potential V-{x) = x^ — 2\x\, 
showed that exact solutions do not exist for the excited states and studied these states by 
a variational method. Finally, starting from the known ground state of V{x) = — 2\x\, 
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we obtained the ground states for the potentials V{x) = and V{x) = x^ + 2|x|, by using 
logarithmic perturbation theory. Comparision with other known results in the literature and 
in the paper are given. 

Some other approaches and improvements can be used to study this class of super- 
potentials. In a forthcoming paper we analyze the solutions for the ground states of 
V±{x) = ± 2|x| by starting with the solutions of V± = ± C and using LPT and 
the 6 expansion of Bender [3| and Cooper |8]. 
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